Relativistic corrections in magnetic systems 
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We present a weak relativistic limit comparison between the Kohn-Sham-Dirac equation and its 
approximate form containing the exchange coupling, which is used in almost all relativistic codes 
of density functional theory. For these two descriptions, an exact expression of the Dirac Green's 
function in terms of the non-relativistic Green's function is first derived and then used to calculate 
the effective Hamiltonian, i.e. Pauli Hamiltonian, and effective velocity in the weak relativistic limit. 
We point out that, besides neglecting orbital magnetism effects, the approximate Kohn-Sham-Dirac 
equation also gives relativistic corrections which differ from those of the exact Kohn-Sham-Dirac 
equation. These differences have quite serious consequences: in particular, the magnetocrystalline 
anisotropy of an uniaxial ferromagnet and the anisotropic magnetoresistance of a cubic ferromagnet 
are found from the approximate Kohn-Sham-Dirac equation to be of order , whereas the cor- 
rect results obtained from the exact Kohn-Sham-Dirac equation are of order 1/c*. We then give a 
qualitative estimate of the order of magnitude of these spurious terms. 
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I. INTRODUCTION 

The relativistic effects play a fundamental role in mag- 
netic systems. Tfiey are responsible for multiple physical 
properties of_great fundamental interest and technologi- 
cal relevanceEl'0: magnetic anisotropy; magnetostriction; 
magneto-optical phenomena such as Faraday effect, Kerr 
effect or magnetic dichroism; anomalous Hall effect and 
anisotropic magnetoresistance (AMR) in metallic mate- 
rials; etc... To describe these properties, one has to use 
either a full relativistic Hamiltonian, i.e. Dirac equa- 
tion, or an effective Hamiltonian, i.e. Pauli equation, 
which includes magnetic interactions and spin-orbit cou- 
pling. However, it is not obvious to find all the different 
terms which have to be taken into account in the effec- 
tive Hamiltonian. For this reason, it appears necessary 
to study the weak-relativistic limit of the Dirac equation 
for magnetic system in order to extract the complete ex- 
pression of the effective Hamiltonian. 

For many-body systems, the relativistic density func- 
tional theory allows to replace the many-body Dirac 
equatior^tuj' the Kohn-Sham-Dirac equation, which has 
the fori 
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where a and /3 are the (4 x 4) Dirac matrices respectively 
related to the (2 x 2) Pauli matrix cr and unit matrix 



a = 



cr 
cr 



(3 = 



1 

-1 
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The effective potential Veff and effective vector potential 
Aeff are functions of external and internal fields as well 
as electron density n{r) = J2t=i and current 

density J(r) = X^ili '0j^(r)caV'i(r): 



n(r') 



r - r' 



dr' 



Sn{r) 



(3) 



Ae^(r) = Aext{r) 



J(r') 



AirefjC J |r 
e 63{r) 



dr' 
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where Vext and Aext are the external potential and vec- 
tor potential and E^c is the exchange-correlation energy 
functional. 

Neglecting orbital cuixents, a Gordon decomposition 
allows to rewrite (|l]) asEfQ 

=c{a:-p)+ (imc^ + Veff - HbP (t • Beff) . (5) 

The last term§l corresponds to the exchange coupling in- 
volving an effective magnetic field Bg^ff = "^ext + ^xc 
where "Qext is the external field and the exchange- 
correlation field given by 
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fJ-B 



i5m(r) 
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where m(r) = X^i^i V'i*'(r)/3cr'0i(r) is the spin density 
and fiB = eh/2m. Expression for the Hamiltonian 
can also be derived and justified with the help of sym- 
metry argumentaJ. 

To describe relativistic effects in magnetic systems, one 
could either start from or from . Even if the sec- 
ond Hamiltonian is an approximate description in com- 
parison to the former one, it has two advantages: on the 
first hand, it is simpler because the vector potential is 
not more present in the kinetic energy and, on the other 
hand, it has a more convenient form for magnetic sys- 
tem because we can include the magnetic interactions 
directly in the scalar potential which is then a sum of 
spin-independent and-, spin-dependent parts. Therefore, 
is generally usedEj t3. However, there is no check on 
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the consistency of these two equations in the weak rela- 
tivistic Hmit. Apart from the orbital magnetism which 
is clearly neglected in , we should obtain the same 
relativistic corrections. We point out in the present pa- 
per that the weak relativistic limits of and dif- 
fer. These differences have quite serious consequences: in 
particular, the magnetocrystalline anisotropy of an uni- 
axial ferromagnet and the AMR of a cubic ferromagnet 
are found from the approximate form to be of or- 
der 1/c^, whereas the correct results obtained from the 
exact form H"^ are of order We then give a qualita- 

tive estimate of the order of magnitude of these spurious 
terms. 

The paper is organized in the following way. In Sec. 
H, we present the derivation of a useful expression of 
the Dirac Green's function in term of the non-relativistic 
Green's function. Then, using this Dirac Green's func- 
tion as a alternative method to separate particles from 
anti-particles, we extract in Sec. Ill the effective Hamil- 



Hamiltonian in terms of (2 x 2) matrices as: 



tonian and effective velocity. These operations are done 
starting both from and in order to compare the 
relativistic corrections of order 1/c^. The calculations 
starting from are given in detail in the text whereas 
the calculations from , which are similar in their prin- 
ciples but lead to different expressions, are summarize in 
the Appendices. Finally in Sec. IV, we discuss the re- 
sults and present a qualitative evaluation of the difference 
between the two descriptions. 



II. DIRAC GREEN'S FUNCTION 

The derivation of the Dirac Green's function has been 
done both for H"^ and . Below, we present in detail 
the calculation for while the calculation for is 
summarize in the Appendix A. In order to simplify the 
notations, we introduce tt = p — eAgff, and write sim- 
ply V instead of Veff. Since we are interested mainly 
in the electrons states in the weak relativistic limit, we 



shift the zero of energy by 



for positrons, we would 



have to shift the zero of energy by +mc . We express the 



V c{(T ■ 7V) 

c{cr ■ Tv) V — 2mc^ 



(7) 



A series of algebraic manipulations allows to express the 
Dirac Green's function G{z) = [z — H^)~^ in terms 
of the (2 X 2) non-relativistic Green's function G{z) = 
[z — H^)^^ associated with the (2 x 2) non-relativistic 
Hamiltonian = {cr ■ iTy/2m + V. We write first G as 
a product of matrices: 



G(z) 



z — c (cr ■ tt) 
-c (cr • tt) 2mc? 

V 

z-V ) ) 

where we have introduced: 

A{z) = 

The explicit inversion of A yields 

A-\z)-- 



Z —c{(T- tt) 

c (cr • tt) 2mc^ 



A-\z), (8) 



(9) 



Goiz) 



■Goiz) 



2mc 

Ti — Gr)(z) Ti — iGoiz) 
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where Gq{z) = (z - (cr ■ ■n:y/2m)-\ Using (|l^ in (|), 
we get 



G(,)=( ^-Go(z)V ^Go{z){z^V) 
^ ' \-§;^Go{z)V l + ^Go{z){z-V) 

^Go{z) §^,Gp{z) 
^ Go (z) Tr^—rGoiz) 



(11) 



The first factor in the right hand side of the above equa- 
tion is calculated by direct matrix inversion, and is equal 
to 



J 



G{z)G^-\z)-G{z)^^{z-V)D{z)^^(g{z)G^ \z)-\) -G{z)§^^{z - V)D{z) 

(g(z)Go^\z) - l) Diz) 

where we have used the relation GqVG = G — Gq and introduced the (2 x 2) matrix 

2mc^ 



D{z)= (l + Q{z)- 



with Q{z) = 1 + (cr ■ 7r)G(z)(cr • Tv)/2m. We report (|T|) in (|Tl|) and finally obtain: 

G(z) ^ ( - ^(^)f^(^ - y)D{z)^A^) G{zmQ-Hz)D{z)Q{z) 

^ ' V m§^A^) ^D{z)Q{z) 



(12) 
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The corresponding expression associated with is 
given in the Appendix A by Eq. ( [A^ ) . We have thus suc- 
ceeded in expressing the (4 x 4) Dirac Green's function in 
terms of (2 x 2) matrices, in particular the non-relativistic 
Green's function G. This formulation provides an alter- 
native method to separate particles form anti-particles in 
the weak-relativistic limit by a simple block diagonaliza- 
tion. More comments and illustration of this method are 



presented in Sec. 

In absence of vector potential (which means tt = p 
in Eq. (|l^ or 6 = in Eq. ([A^)), the method and re- 
sult, are similar to those presented earliep-bji, Gesztesy et 
alxA and more recently by Brouder et a/.HIij for the non 
magnetic case, although the form we give here is simpler 
than the one given by these authors. The generalization 



to magnetic cases suggested in Refs. |1J,|19| is not correct 
because the authors, who have started from , have 
omitted the factor (3 in the expression of the exchange 
coupling -^B/3(cr • Bcff). 

In the limit of low electron energies, it is useful to have 
a semi-relativistic expression. In order to get it, we per- 
form an expansion of the operator D in powers of 1/c 
which allows to write the Dirac Green's function as a 



series 



G{z) 



(15) 



where G^"'\z) is the term of order 1/c". The successive 
terms are 



G(°)(z) = 
G(2'=+i)(z) = 

G(2"+2)(z) = 



G{z) 




G(z) 



2m c V 2mc^ 





'G{z)q;^A-~V)(~Q{zY-^) §^G(z) 



(z-V) 
2mc^ 



Q{z) 



(16) 
(17) 

(18) 



for fc > 0. We remark that odd terms in the expansion 
of the Green's function in powers of 1/c are odd matrices 
whereas even terms are even matrices. Corresponding 
expressions for the Dirac Green's function expansion as- 
sociated with are given in Appendix A. 

One advantage of the expression ( [Tsl ) is that we can di- 
rectly identify and calculate the terms which gives rise to 
particular effect according to the order of this effect with 
1/c. Consider for example the magnetic anisotropics. 
In the case of a system with uniaxial anisotropy, the 
anisotropy energy is quadratic in the spin-orbit coupling 
\so (i.e., of order 1/c^). As the anisotropy energy is linear 
with respect to the Green's function, we have to consider 
the term G^'^^ For cubic anisotropy, the anisotropy en- 
ergy is of order Af^ (i.e., of order 1/c^), then we have 
to consider the term G^^\ For galvanomagnetic and 
magneto-optical effects, one needs to calculate the con- 
ductivity tensor. The latter is expressed as a product of 
two Green's functions and two velocity operators, which 
are cot in the relativistic theory. For the effects which are 
linear in spin-orbit coupling, (i.e., of order 1/c^), such 
as the anomalous Hall effect, and the Kerr and Faraday 
magneto-optical effects, one needs to calculate all terms 
up to G'^^'^ . For the effects that are quadratic in spin-orbit 
coupling (i.e., of order 1/c'*), such as the AMR or the 
magnetic birefringence, one needs all terms up to G^^^ 



Note that the usual (i.e., non-relativistic) conductivity is 
obtained from the terms up to G^^^ . 



III. WEAK RELATIVISTIC LIMIT 

A. Effective Hamiltonian 

We turn now our attention to the effective Hamil- 
tonian which can be obtained through different ways 
such as for example the elimination of the lower com- 
ponents of the wave function or the Foldy-Wouthuysen 
transformation |Wliich requires a succession of canonical 
transformationsE2l. As we have the explicit expression 
of the Dirac Green's function in term of (2 x 2) matri- 
ces, we do not need to used such methods to separate 
the particles from the antiparticles. Actually, a block di- 
agonalization of (^^ allows to cancels the terms which 
couple the upper and lower components and then to ex- 
tract the effective Hamiltonian. Let us start with the 
non-relativistic limit (c — > oo) which is obtained in a 
transparent manner. Indeed, in this limit, the matrix el- 
ements of G = G'"' are different from zero only in the left 
upper part (see Eq. (14)). Then the separation between 
particles and antiparticles is naturally made: the Green's 
function which described the particle are directly done by 
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G and the effective Hamiltonian is H"^. In the general 
case (arbitrary value of c), the separation between parti- 
cles and antiparticles can be exactly made only for free 
electrons. For particles in a potential, like in our case, a 
development in powers of 1/c has to be performed which 
means a restriction to the weak-relativistic limit. Indeed, 
the block diagonalization of the Green's function is only 
possible if we cut the expansion (|l^) . S-uch limitation is a 
common characteristic for all methods^. Since we restrict 
our calculations to the lowest order relativistic correc- 
tions to the Hamiltonian, we have to cut the expansion 
(H) of the Gr een's function after the second order with 
1/c: 



diagonalization of H"^). This result is exactly similar to 
that given in the literature in the case of a single parti- 
cle Dirac equation (see for example Ref. In order to 
get a more usual expression of i?^, we have to perform 
some transformations which are detailed in Appendix B. 
In the case of an uniform effective magnetic field and ne- 
glecting orbital magnetism, H^^ reduces to (see Eq. ( p^ ) 
in Appendix B): 



P 



AT/- 



iyv X tt) 



(25) 



G(z) 



2mc y f 2mc'^ ^ V / 



(19) 



The block diagonalization of G corresponds to the change 
of basis M^^GM where the (4 x 4) unitary matrix M is 
given by: 



M 



1 + ^ 



2mc 

2mc 



and leads to the block diagonal Green's function: 



M-^Giz)M ^ { 



(20) 



(21) 



By means of this transformation, wc have achieved a de- 
coupling between the particles and the antiparticules: the 
(2 X 2) matrix describes the particles whereas the 
(2 X 2) matrix g_ describes the antiparticles. We get: 



- cr • TV T 



G{z 



Sm'^c'^ 



and 



2mc 



2mc 



''g(z), 



(22) 



(23) 



Since now, we restrict our study to the particles. Using 
the relations Gz = 1-1- GH^ = 1-1- H^G, we can trans- 
form the last term in (^2|) and write under the form 
G + GH±G where: 



We obtain the usual relativistic corrections (relativistic 
mass correction, Darwin term and spin-orbit coupling) 
plus an additional contribution due to the presence of the 
exchange coupling. For a non-uniform effective magnetic 
field, which is the case in realistic problems, further rel- 
ativistic corrections are obtained (see Eq. (B4)). Similar 
calculations, presented in Appendix C, have been done 
starting from by performing a block diagonalization 
of the Dirac Green's function calculated in Appendix A. 
We can now compare the relativistic corrections iJ^ and 
H^^ obtained in the two descriptions in the case of a uni- 
form effective magnetic field. Comparing ( p5| ) and (|C^), 
we observe two differences: one in the spin-orbit coupling 
because the non-relativistic velocity v is equal to tt/to in 
the first description whereas it is equal to p/m in the 
second description; and an other one in the relativistic 
corrections Hrxc to the exchange coupling, given by the 
last term in Eqs. (25) and (C7), respectively. Actually, 
from d23) we have: 



whereas from ( |C7| ) we have 



P ('^•Beif), 



(26) 



(27) 



What is problematic is that iJ^^ and H^^^ couple the 
spin and momentum in a quite different manner. Com- 
ments and consequences of this difference are presented 
in the section IV. 



B. Effective velocity operator 



(^. 



1 



(^. 



TT 



-V -V 



rj(cr • TT)V{a- ■ tt) 



(24) 



This expression corresponds to the relativistic corrections 
of order 1/c^ to the non-relativistic Hamiltonian H^. 
Thus, the effective Hamiltonian, i.e. Pauli Hamiltonian, 



is H^ff 



+ H^^ (it can also be derived from a block 



To complete this study, we want to comment briefly 
on the velocity operator which is by definition: 

1 1 f)]-f 

v=-[r,H] = -{rH-Hr) = —, (28) 
in in ap 

where r is the position operator. When we report the 
expression of the Hamiltonians or ^ we get the 
simple form: 
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co- 
co- 



(29) 



The effective velocity can be obtained from Eq. (2^ by 
the change of basis M~^vM but in order to get the cor- 
rections of order 1 /c^ to the velocity, it would be neces- 
sary to expand M up to the order 1/c^ which is cumber- 
some. It can also be obtained form H^jj using : 



'eff 



ih 



effi 



in in 



(30) 



where v'^ is the non-relativistic velocity and v^^ the rel- 
ativistic corrections of order 1/c^ to the velocity. From 
the expression of H^^ we get = ir/m and from (B4) 
where diamagnetic terms are neglected, we get: 



MB 

2m^c^ 

e 



:(o- X VV) 



(P('^-Beif) + (o--Bejy)p) 

(p^A + Ap^ + p(p • A + A • p) 



+ (p- A + A-p)p). 



(31) 



In the case of a uniform effective magnetic field and in 
absence of orbital magnetism, it reduces to: 



fJ-B 



(o- X VV) 



p(o- • Beff ) 



(32) 



The first term is the contribution which comes from the 
relativistic mass correction. The second term, the so- 
called anomalous velocity, results from the spin-orbit cou- 
pling and can play an important role, for example, it 
leads to the side-jump mechanism in the anomalous Hall 
effect. The last term is due to the presence of the ex- 
change coupling and has no specific name. Its symme- 
try is also different in comparison to the relativistic cor- 
rections to the velocity obtained starting from (see 
Eq. (PI) in Appendix D). 



IV. DISCUSSION 

This study gives some clarification concerning the as- 
sumptions made when one replaces the Hamiltonian 
by the Hamiltonian . Even if the consequences of 
such approximation are not fully known as Kiibler justly 
noticesEj; it was generally believed that the transforma- 
tion from to neglects only orbital magnetism 
effects. However, the calculations of the weak-relativistic 
limits of H"^ and made in Sec. HI, reveal an ad- 
ditional difference which corresponds to a different sym- 
metry of the relativistic corrections to the exchange cou- 
pling: whereas -ff^cc (see Eq. (p6|)) is isotropic with 



respect to the direction of the momentum, -ff^^, (see 
Eq. (p?!)) is anisotropic because its amplitude depends, 
through the scalar product (p-Bgjff ), of the angle between 
the directions of the momentum and the effective mag- 
netic field. Thus, the use of can lead to anisotropic 
effects, such as the magnetocrystalline anisotropy or the 
AMR, which differ from those obtained from the exact 
Kohn-Sham-Dirac Hamiltonian H^. 

Let us first consider the magnetocrystalline anisotropy 
of a uniaxial system (e.g., a material with an hexagonal 
lattice, or an ultrathin film). If we start from H^, the 
magnetic anisotropy arises only as a second order per- 
turbation due to the spin-orbit coupling, so that it is of 
order 1/c^. In contrast, if we use , it is easy to see 
that the relativistic correction of the exchange interac- 
tion, gives rise to an additional contribution to the 
magnetocrystalline anisotropy, already in the first order 
of perturbation, i.e., of order 1/c^, which is unphysical. 

Although the details are somewhat more complicated, 
a similar result is obtained when considering the AMR of 
a cubic system: starting from H^, the AMR arises as a 
second order perturbation due to the spin-orbit coupling, 
i.e., it is of order 1/c'*, whereas starting from H^, an ad- 
ditional AMR term of order 1/c^ arises as a first order 
perturbation due to 



rxc 




10" 10'" 10' 
c (m.s"^) 



10'- 



FIG. 1. This figure shows the typical variation of the 
magnetocrystalline anisotropy of a uniaxial ferromagnet as a 
function of c. The physical value of c is indicated by the ver- 
tical line. The solid curve is the correct value obtained 
from H"^, whereas the dashed line represents the spurious 
contribution K obtained from H . 



In the limit of c ^ oo, these spurious terms of order 
1 /c^ dominate over the correct terms of order which 
is of course unacceptable. As hII, ah initio calculations of 
magnetocrystalline anisotropyOliS starting for the Dirac 
equation rely on the approximation (^, there validity 
can be questioned a priori. Let us make a simple esti- 
mate of the orders of magnitude for the physical value of 
c.|-Starting from H^, the magnetocrystalline anisotropy 
isE^ 
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W 



(33) 



where Ago is the spin-orbit constant and W the band- 
width. Starting from one other hand, it is quite 
easy to calculate the magnetocrystalline anisotropy due 
to M 



rxc 



K 



B 



where is the Fermi level and Aea^ the exchange split- 
ting. For a transition metal ferromagnet, by taking typ- 
ical values Xso ~ 0.1 eV and W ^ 5 eV, one obtains 
K 2x 10-3 eV. Taking ep « 10 eV, Ae^ « 2 eV, 
and mc^ « 500 keV, we obtain ^ 2 x 10"^ eV. 
Therefore, in spite of the fact that ^ in the 
limit c ^ cxD, we find that <^ for the the 

physical value c « 3 x 10^ m.s^^. This result is vi- 
sualized in Fig. 1. Therefore the quantitative results 
of first-principles calculations of the magnetocrystalline 
anisotropy based on the approximate Kohn-Sham-Dirac 
Hamiltonian should not be perturbed in a signifi- 
cant manner by the spurious contribution of order 1/c^. 
In spite of this fortunate circumstance, it would be desir- 
able to develop a more satisfactory theoretical approach 
which is free from unphysical spurious contributions. 



For a non-uniform effective magnetic field, additional 
relativistic corrections are obtained (see Eqs. (B4) and 
(|C6D). The qualitative importance of these corrections 
remains to be investigated. 

To conclude, we want to underline two aspects which 
represent the important results of this work. On the one 
hand, we have derived a convenient form of the Dirac 
Green's function which is valid for any value of c and has 
the advantage to be express in term of (2 x 2) matrices, 
in particular the non-relativistic Green function. This is 
a quite general form which can be used to treat differ- 
ent problems. We have applied it as an alternative way 
to separate particles from anti-particles in order to ex- 
tract the effective Hamiltonian and effective velocity in 
the weak-relativistic limit. On the other hand, we have 
performed a detailed comparison in the weak-relativistic 
limit between the Kohn-Sham-Dirac equation and its ap- 
proximate form containing the exchange coupling. This 
study has revealed a difference of symmetry in the rel- 
ativistic corrections to the exchange coupling which can 
leads to artificial anisotropic effects. A qualitative esti- 
mation has shown that, in the case of transition metals, 
this difference is not significant. 



APPENDIX A: DIRAC GREEN'S FUNCTION IN 
PRESENCE OF AN EXCHANGE COUPLING 



The Dirac Green's function G(z) 
be written as: 



(z — H^) -'^ can then 



In this appendix, we summarize the derivation of the 
Dirac Green's function starting from given by (||). 
It follows the same steps than starting from but in- 
volves different matrices and leads to a different expres- 
sion of the final Dirac Green's function. To simplify the 
notations, we introduce O — ~fiB ("" • '^eff), then: 



V + 
c (cr 



e 
p) 



V 



C((T • p) 

— 2jnc? - 



e 



(Al) 



G{z) = [l + A-\z) 
xA-\z). 



-V -Q 
z~V+Q 



(A2) 



is given by (|T^) when we replace tt by p and define 
Gq{z) = {z — / 2m)~^ . We perform the direct inversion 
of the matrices which appear in ( AS ). It leads to the final 
expression: 



G(z) = 



( G{z) - G{z)^(z -V + Q)D{z)^G{z) G{z)^0-\z)D{z)Q{, 



2mc' 



h^D{z)Q{z) 



(A3) 



where G is the (2 x 2) non-relativistic Green's function 
associated with the (2 x 2) non-relativistic Hamiltonian 
j^B _ p'i I2jn -JrV + Q, and the operators D and Q are 
given by: 



D{z)= 1 + Q(z) 



{z-V + Q) 
2mc^ 



Q{z) = 1 + 



(<T.p)G(z)(a.p) 

2m 



(A5) 



A semi-relativistic expansion G{z) = X]^o (^) '"^^ 
(^^) also be given. The successive terms are 
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for fc > 0. 



G(z) 






2mc^ j 2n 



:G(z) 



-Giz)§^Az-V + e){-Q{z)^^^) ^G{z) 



(A6) 
(A7) 

(A8) 



APPENDIX B: EFFECTIVE HAMILTONIAN IN 
PRESENCE OF A POTENTIAL VECTOR 

Transformations of and H:^^ allow to get a more 
usual expression of the effective Hamiltonian H'^g. Ne- 
glecting diamagnctic terms, we have: 

-^(p^A + A^p). (Bl) 

where we have used the identities: 

(cr • tt) (cr • tt) = TT • TT + zcr • (tt X tt) , 

pxA + Axp = —iKBeff- (B2) 

For a uniform effective field, reduces to: 



(B3) 



where L = r x p is the orbital momentum. In absence 
of orbital magnetism, this last expression is identical to 
H^. We turn now our attention to the relativistic cor- 
rections -ff^. After transformation of (|2j) and neglecting 
the diamagnetic terms, we can rewrite H^^ as: 



P 



AV- 



a ■ iyV X tt) 



+ 



/ 2 2 

+ ( — (P •A + A-p) + (p-A + A- p)^ 

+ (A.p-p.A)y-V^(A.p-p.A)), (B4) 
where we have used the following identities: 

{(T ■ tt) V {cr ■ tt) = ttV ■ tv + icr ■ (nV x tt) , 

2pV ■p=p'^V + Vp^ + h^AV, (B5) 

and replaced the momentum operator p as well as the 
operator p'^, when they act only on the potential, respec- 
tively by the gradient V = —p/ih and by the Laplacian 
A = ~p^/h^. In the case of a uniform effective magnetic 
field, H^^ reduces to: 



AV - 



Arrive 



J(T ■ {W X tt) 



pHi<T + L)-B,ff). 



(B6) 



APPENDIX C: EFFECTIVE HAMILTONIAN IN 
PRESENCE OF AN EXCHANGE COUPLING 

In this appendix, we summarize the derivation of the 
effective Hamiltonian starting from the Dirac Green's 
function (A3). We cut the expansion after the second 
order with 1 /c, then: 
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The block diagonalization M ^GM is obtained with: 

1 



A/ = 1 
and leads to: 

g+{z) = G{z) + 



p 



TP 

2mc 



P 

2nic 
1 



(C2) 



:G{z) + G{z) 



'G{z)^-^{z -V + e)^G(z), 
^ ' 2mc ^ ' 2mc ^ ' 



(C3) 



which can be written under the form G + GH^^G where: 

.4 1 



H 



B 



P 



^2 



• p) (y - 9) ((T • p) 



p 



(v + e)- {v + e) 



p 



(C4) 



In order to get a more us ual expression of , w e made 
some transformations of (C4) using identities (|B5|) where 
we have replaced tt by p and the relation: 

(<T-p) (cr-Be#)(cr-p) = 

- p (O- ■ Beff) • p + a (V X Beff) ■ P 

+ i(T-p) {B,ff • p) + (p • B,^) (cr • p) 
Finally, we get: 

"4 



(C5) 



rrB _ P 
J^^.. — — 



, , A{V-^iB{(T-B,ff) 
-^^■iVV.p) + ^iV><B.,).p 



fJ-B 



(('^•p)(Beif p) + (p-B,^)(^.p; 



(C6) 



This expression of the effective Hamiltonian in presence 
of the exchangfi^oupling is equivalent to the one obtained 
by Kraft et aZo. In the case of a uniform effective mag- 
netic field, it reduces to: 



P 



+ -^a ■ {VV X p) 
(<T.p)(B,^.p). (C7) 



APPENDIX D: EFFECTIVE VELOCITY 
OPERATOR IN PRESENCE OF AN EXCHANGE 
COUPLING 

In this appendix, we summarize the derivation of the 
effective velocity starting from . From the expression 
of = /2m + y + 0, we get v"^ = p/m and from 
the expression ( |C^ ) of i?^, we get: 



+ i£|2HBeifP) + (p-B.^) 

+ {^-p)B,ff+B,ff{cT-p)). (Dl) 



This expression differs from the result obtained by Kraft 
et alx^ on two points: they get a wrong coefficient (a 
factor 2 missing) for the contribution of the relativistic 
mass correction and they obtain a contribution from the 
Darwin term which should not appear because the com- 
mutator [r, AT/] is equal to zero. In the case of a uniform 
effective magnetic field, (Dl) reduces to: 



rc 



(D2) 
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